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A systematic investigation of the average multiplicities of light charged particles and intermediate 
mass fragments emitted in peripheral and semiperipheral collisions is presented as a function of 
the beam energy, violence of the collision and mass of the system. The data have been collected 
with the FIASCO setup in the reactions 93 Nb+ 93 Nb at 17, 23, 30, 38AMeV and 116 Sn+ 116 Sn at 
30, 38 A MeV. The midvelocity emission has been separated from the emission of the projectile-like 
fragment. This last component appears to be compatible with an evaporation from an equilibrated 
source at normal density, as described by the statistical code Gemini at the appropriate excitation 
energy. On the contrary, the midvelocity emission presents remarkable differences for what concerns 
both the dependence of the multiplicities on the energy deposited in the midvelocity region and the 
isotopic composition of the emitted light charged particles. 

PACS numbers: 25.70.-z, 25.70.Lm, 25.70.Pq 



I. INTRODUCTION 

In recent years many works (see, e.g., 0, 0, 0, S 
and references therein) were devoted to the investigation 
of Light Charge Particle (LCP) and Intermediate Mass 
Fragment (IMF) emission in semiperipheral and periph- 
eral heavy ion collisions at Fermi energies. It is well 
established that these reactions show mainly a binary 
character, with two heavy remnants in the exit channel, 
possibly undergoing a subsequent sequential fission |(|. 
Among the various features of the emission, one in par- 
ticular raised a lot of interest: a large amount of emission 
is located at midvelocity (i.e., close to the center-of-mass 
velocity), mainly for the IMFs but also for lighter parti- 
cles [1 U H II IE 03 The origin of this phenomenon 
is still actively debated. Possible interpretations range 
from, e.g., a kind of neck rupture 0, ^J, due to me- 
chanical and/or chemical instabilities, to a fast statistical 
emission from one of the heavy fragments perturbed by 
the proximity of the other heavy remnant |l2l If3 | . 

In a previous work by means of a three-body 
Coulomb trajectory calculations, we put into evidence 
two components in the experimental IMF emission pat- 
tern: a fast emission from the phase space region in be- 
tween the two heavy remnants, somewhat reminiscent of 
a neck fragmentation or a participant-spectator model, 
and a later emission from the (possibly deformed) sur- 
face of the heavy fragments. The latter mechanism may 
be interpreted as the evolution of the fast oriented fission 
mechanism 0, IbH Il5| toward very large mass asymme- 
tries. Similar conclusions were drawn in a more recent 
paper |l6j |. using the same basic approach of Coulomb 
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trajectory calculations. 

The midvelocity emission is particularly evident in 
very peripheral reactions (b/b graz > 0.9), where the kine- 
matics of the collision is relatively simple and the mul- 
tiplicity of LCPs (and IMFs) emitted from the two fully 
accelerated main fragments is quite small (about 0.3 LCP 
per event in Ref. |3j). As a consequence, the investiga- 
tion of peripheral collisions is a fundamental tool in order 
to shed more light on the production mechanism and on 
its evolution with the beam energy. The short interac- 
tion times typical of these collisions can also be used to 
investigate isospin diffusion and equilibration in collisions 
between nuclei with different N/Z ratio (see an d ref- 
erences therein). 

In Q we found that in peripheral collisions the en- 
ergy stored inside the midvelocity matter and the ex- 
citation energy of the quasi-projectile have similar val- 
ues. This evidence suggested a larger energy density (> 7 
MeV/nucleon) in the midvelocity "source", with respect 
to the excited quasi-projectile (< 2 MeV/nucleon), since 
the mass localized at midvelocity is certainly smaller (the 
value would even rise to «13 MeV/nucleon if the mass 
of the midvelocity "source" were identified with the to- 
tal mass of its emitted particles). Similar conclusions 
were drawn from the comparison of transverse velocities 
too 01- If one assumes, as stated for example in [l^ . 
that in central collision the multifragmentation process 
starts at excitation energies of the source greater than 
«3AMeV, then the midvelocity emissions might be in- 
terpreted as a first appearance of the multifragmentation 
process ^lj. Indeed it was also claimed that midveloc- 
ity fragments associated with midperipheral and central 
collisions present similar characteristics |l8j . 

This paper presents a systematic investigation of the 
average multiplicities of LCPs and IMFs emitted in pe- 
ripheral and semiperipheral collisions of symmetric sys- 
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terns at Fermi energies. The evolution of the multiplic- 
ities is studied as a function of the excitation energy of 
the emitting "source" , mass of the system and beam en- 
ergy. Section [H] briefly describes the experimental setup 
FIASCO [2(j and the investigated reactions. Section ITTll 
discusses the event selection and describes the analysis 
methods for separating the evaporative and midvelocity 
components of the emissions. Scction llVl presents the ob- 
tained results on the particle multiplicities and their scal- 
ing with the excitation energy and mass of the "source" 
and the beam energy while conclusions are drawn in Sec- 
tion |V| Some more technical points are presented in the 
Appendices. 



TABLE I: Some calculated reaction parameters 

93 Nb+ 93 Nb 116 Sn+ 116 Sn 
E/A (MeV) 17 23 30 38 30 38 



E ; c n m - (MeV) 


791 


1069 


1374 


1772 


1715 


2210 


v rc i (mm/ns) 


57.3 


66.6 


75.6 


85.8 


75.6 


85.8 


Olrt (deg) 


7.8 


5.5 


4.2 


3.2 


4.8 


3.6 


fcgraz (fm) 


11.2 


11.6 


11.9 


12.1 


12.5 


12.8 


^graz (Jl) 


470 


568 


660 


762 


865 


1002 


calc / i \ 
O-rcac (mb) 


3938 


4260 


4462 


4621 


4942 


5143 



II. EXPERIMENTAL SETUP 

The results presented in this paper were obtained in 
a systematic study 0, 0, 0, ^ of heavy ion collisions 
performed at the Superconducting Cyclotron of the Lab- 
oratori Nazionali del Sud of INFN in Catania. 

Beams of 93 Nb at 17.0, 23.0, 29.5 and 38.1 A MeV and 
of 116 Sn at 29.6 and 38.1 AMeV were used to bombard 
93 Nb and 116 Sn targets of about 200 ^g/cm 2 thickness. 
This paper concerns solely the symmetric collisions 93 Nb 
_l_ 93 Nb and 116 Sn + 116 Sn, for which some parameters, 
calculated using the Bass interaction distance Rj n t [23[ . 
are listed in Table H 

The experimental data have been measured with the 
FIASCO setup (described with more details in [2(1 l2ll|). 
a multidetector particularly well suited for the study of 
peripheral and semiperipheral reactions, where only few 
heavy remnants are produced. In fact, as a character- 
istic feature, this setup includes 24 large-area position- 
sensitive Parallel Plate Avalanche Detectors (PPAD) 
0, 13 covering about 70% of the forward hemisphere. 
They are fully efficient for heavy fragments with Z>10 
and are used to measure velocity vectors with very low 
detection thresholds (~ 0.1 AMeV) and good accuracy 
(position and time-of-flight resolutions of 2-4 mm and 
700 ps (FWHM) , respectively, with flight-paths of about 
3.5 m for 9 < 10°). In this way it is possible to simul- 
taneously detect both the projectile-like fragment (PLF) 
and the very slow target-like fragment (TLF) even in the 
most peripheral collisions and to perform a kinematic re- 
construction of the events j2f| . 

Behind the most forward PPADs, a mosaic of 96 Silicon 
telescopes allows to measure the energy, the charge, and 
the final mass (by means of the time-of-flight) of the PLF. 
Each telescope (with 28 x 28 cm 2 active area) consists of 
a 200/im AE detector followed by a 500^tm E res detector. 

Finally, the setup is completed by 182 three-layer 
phoswich scintillation telescopes, mounted behind most 
of the PPADs and covering about 30% of the forward 
hemisphere (plus a reduced sampling in the backward 
hemisphere), which are devoted to the detection of light 
charged particles and intermediate mass fragments. The 
phoswiches allow isotopic identification of Hydrogen (and 



in some cases of Helium) isotopes and charge identifica- 
tion of heavier products up to Z ^15-20 (with a thresh- 
old of about 3-3.5 A MeV) and give a direct measurement 
of the time of flight (and hence of the velocity) of all these 
reaction products, without any need for tricky and time- 
consuming energy calibrations of the scintillators. 

The experiment was performed together with the ho- 
doscope Hodo-Ct [2(J of the Temperature experiment 
|27j . but its data have not been used here. 

The results presented in this paper are focused on bi- 
nary events - by far prevailing in peripheral and semipe- 
ripheral collisions - with only two large reaction remnants 
(Z>10) detected by the gas counters, and on the associ- 
ated multiplicities of LCPs and light IMFs with Z=3-7. 

III. EXPERIMENTAL DATA 
A. The choice of an "ordering parameter" 

For a meaningful sorting of the data in homogeneous 
bins of increasing centrality, one needs an experimental 
variable which is expected to have a monotonic and pos- 
sibly narrow relationship with the impact parameter of 
the collision. This "sorting" variable should also be -as 
much as possible- independent of the studied observables, 
in order to avoid (or reduce) autocorrelation effects. 

Global variables (like multiplicities, flow angles, trans- 
verse energies, etc.) are not used in the present work, as 
they are best suited for experiments covering very large 
solid angles (close to Air). Therefore other variables, re- 
lated to more particular aspects of the reaction, must be 
considered. Since the main subject of this paper is the 
emission of LCPs and IMFs in binary collisions, it is nat- 
ural to restrict the choice to variables which make use 
of experimental information concerning the two main re- 
action partners (like, e.g., the secondary charge Z scc of 
the PLF residue, or its secondary velocity «plf, or the 
relative velocity u re i between PLF and TLF). 

Our choice is a variable [2^ defined as 

TKEL = E; c n m ' — lA^reb (1) 

where E; c n m - is the center-of-mass (cm.) energy of the 
collision in the entrance channel, i> re i the reconstructed 
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relative velocity and fl the reduced mass calculated with 
the masses obtained from the kinematic analysis. 

It is to be noted that -by definition- the kinematic 
method constrains the primary masses of the two reaction 
partners to add up to the total mass of the system. Thus, 
while at low incident energies, where reactions are strictly 
binary, TKEL may truly represent the "Total Kinetic En- 
ergy Loss" of the collision - namely the total amount of 
kinetic energy transferred from the relative motion into 
internal energy of the colliding system - it is important 
to note that at Fermi energies, where the presence of a 
sizable midvclocity emission causes an overestimation of 
the kinematically determined ft, this interpretation is no 
longer correct. Therefore, in this work TKEL is used just 
as an ordering parameter for sorting the events in bins 
of increasing centrality. 

A short discussion about the relationship of TKEL 
with other possible sorting variables is also presented 
in Appendix A, together with an estimation of the cor- 
respondence between TKEL and impact parameter ob- 
tained both by means of model calculations and exper- 
imentally via a direct integration of the reaction cross 
section. We anticipate here that, for a given colliding 
system, equal values of TKEL correspond to similar im- 
pact parameters regardless of the bombarding energy, at 
least as far as (semi-)peripheral collisions are concerned. 
When the system is changed (in our case from Nb+Nb 
to Sn+Sn), a given value of TKEL indicates a somewhat 
larger impact parameter for the heavier system, as ex- 
pected from the larger nuclear radii. 



(more details arc given in Appendix A). 

The third row shows that the exit channel with two 
heavy remnants remains the dominant one even at the 
highest investigated energies, where it still accounts for 
more than half of cr^^, thus demonstrating the persis- 
tence of the binary or quasi-binary character of the re- 
actions. The three-body channel, which is likely to be 
due to sequential fission or to the "fast oriented fission" 
^3 °f onc °f the two main fragments, adds an ap- 
preciable contribution. For the 93 Nb beams (where more 
beam energies have been studied) it is interesting to note 
that at 17 A Me V the two- and three-body channels ac- 
count for nearly 100% (within errors) of the whole reac- 
tion cross section and that this percentage decreases with 
increasing bombarding energy. This behavior may be 
due to the progressive opening of new reaction channels 
(four- or more-body reactions, possibly multifragmcnta- 
tion). However, even at the highest energy of 38AMeV, 
the two- and three-body exit channels altogether still rep- 
resent about 75% of 

The main uncertainty on the quoted numbers concerns 
the cross sections for two-body events. It mainly arises 
from the difficulty of separating clastic and quasi-clastic 
scattering in the vicinity of the grazing angle. The ex- 
perimental data have been integrated starting from the 
calculated grazing angles of Table HJ after having verified 
that they agree with the experimental "quarter-point" 
angles within a few tenths of degree. Such an uncer- 
tainty in angle corresponds to an uncertainty of the order 
of 200-300 mb on the cross sections and of about 5-6% 
on the percentages in the last row. 



B. The cross sections 

For each reaction, the simultaneous measurement (via 
a dedicated and suitably down-scaled "singles" trigger) 
of elastically scattered projectiles hitting the most for- 
ward gas detectors has been used for determining the 
conversion factor ( "millibarn-per- count" ). After correct- 
ing for the experimental filter, this allows to perform a 
quantitative estimate of the experimental cross sections 
pertaining to the detected two- and three-body events. 
They are summarized, for all investigated reactions, in 
the first and second row of Table [HI while the third row 
indicates the percentage of the calculated total reaction 
cross section which is accounted for by these two channels 



TABLE II: Reaction cross sections 







93 Nb - 


- 93 Nb 




116 Sn+ 


116 Sn 


E/A (MeV) 


17 


23 


30 


38 


30 


38 


a 2b P ( mb ) 

°3b P ( mb ) 
exp / calc 
2b / rcac 
cxp / calc 

^b+ab/^roac 


3210 
640 

82% 
98% 


2840 
940 
67% 
89% 


2750 
940 
62% 
83% 


2810 
790 
61% 
78% 


2670 
1030 

54% 
75% 


3010 
980 

58% 
77% 



C. The emission pattern 

The shape of the emission pattern of LCPs and IMFs 
can be best appreciated from the distribution of the ex- 
perimental yields in the (v»,v±) plane, where v\\ and 
v± are the velocity components (in the cm. refer- 
ence system) parallel and perpendicular, respectively, to 
the asymptotic PLF-TLF separation axis (for the TKEL 
range addressed in this work, the separation axis remains 
in most cases within about 10° from the beam axis). As 
the solid angle coverage of the FIASCO setup, although 
large, is however significantly smaller than 47r, it is first 
necessary to correct the data for the limited geometrical 
coverage of the setup and for the low-energy identifica- 
tion thresholds |2ll| 

The correction 0, [22] is obtained from a Monte Carlo 
simulation which produces - for each particle species - a 
random isotropic emission. In fact, it has been verified 
that, thanks to the large acceptance and axial symmetry 
of the setup, the obtained correction is largely indepen- 
dent of the specific emission pattern used in the sim- 
ulation. Therefore the emission adopted in the simula- 
tion homogeneously fills a sphere in phase space, centered 
in the cm. origin with radius t | ma x = 120 mm/ns, thus 
covering all regions which are relevant for the processes 



4 



TKEL=200MeV TKEL=500MeV TKEL=800MeV TKEL=200MeV TKEL=500MeV 




V, (mm/ns) 

E/A=38MeV E/A=23MeV 11 

FIG. 1: (color) Experimental yields (corrected for the setup efficiency) of protons (upper row), a particles (middle row) and 
IMFs with Z=3-7 (lower row) for 93 Nb + 93 Nb at 38 (left panels) and 23AMeV (right panels). Yields - with logarithmic 
level spacing in factor-of-2 steps - are in the plane (vn,v±), with respect to the PLF-TLF separation axis, the origin being 
in the center-of-mass system. Columns refer to 100-MeV wide TKEL bins centered, from left to right, at 200, 500, 800 MeV 
(at 38 A MeV) and 200, 500 MeV (23 A MeV). The circles in the TKEL=200 MeVpanels show regions affected by velocity 
thresholds (due to the thin fast-plastic scintillators of the phoswich telescopes, see |2fij[1. 



under study. The velocity vector of each LCP or IMF 
is then transformed into the laboratory reference frame 
and the appropriate experimental filter of the FIASCO 
setup is applied (geometry and identification thresholds 
of the phoswich detectors). In order to mimic as closely 
as possible the binary reaction step, all its relevant pa- 
rameters (like the orientation of the PLF-TLF separation 
axis and the PLF and TLF cm. velocities) are those ob- 
tained with the kinematic reconstruction directly from 
the measured two-body events. Using these parameters, 
the procedure allows to determine, for each cell in the 
(vu, v±, </>) space and for successive bins in TKEL, an 
average correction factor, which is obtained as the ra- 
tio between the number of generated particles and the 
number of particles surviving the experimental filter of 
the setup. (Here <j> is the out-of-plane angle with respect 
to the reaction plane; this coordinate was explicitly con- 
sidered for a better efficiency correction in case of pos- 
sible out-of-plane anisotropies in the experimental data, 
e.g., due to angular momentum effects). When analyzing 
the actual data, the experimental yields of each particle 
species are multiplied, cell-by-cell, by the now described 
correction factors. 

As an example of the obtained distributions, Fig. ^ 
shows the efficiency-corrected experimental yields of pro- 
tons, a particles 1 and IMFs with Z=Z-7 (first, second 



1 Among the Helium isotopes, 4 He is the dominant one, while 3 He 
accounts for a few percent of the total (it is visible on the left tail 
of 4 He only in the phoswiches with the best resolution) and 6 Hc 
cannot be singled out. Therefore, although the data presented 
in this paper refer to all He isotopes, they are representative of 
the behavior of a particles only. 



and third row, respectively) in the reaction 93 Nb + 93 Nb 
at 38 and 23 A MeV for three and two bins of TKEL, 
respectively. In this presentation, in absence of instru- 
mental effects, the average positions of the PLF and 
TLF emitters lie -by definition- on the horizontal axis, 
symmetrically with respect to the cm. origin. One 
clearly sees characteristic circles around the positions of 
PLF and TLF, indicating a Coulomb-dominated emission 
from these sources, as it is expected for the sequential de- 
cay of such highly excited systems. (In the case of TLF, 
the inner part of the Coulomb circles is marginally af- 
fected by the velocity thresholds.) It has also to be noted 
that, since yields are presented in Fig.^- and not invari- 
ant cross sections - the intensity of an isotropic emission 
must gradually decrease along the Coulomb circles while 
approaching the v ii-axis, until it vanishes when it reaches 
this axis. At parallel velocities intermediate between 
those of PLF and TLF, one observes an additional contri- 
bution, the so-called "midvelocity" (or "neck") emission. 
Although present for all particle species, this emission is 
particularly evident and important for a particles and 
even more for IMFs. 

The FIASCO setup has a much better solid angle cov- 
erage in forward direction, where -in addition- thresholds 
effects do not play any practical role, because the ener- 
gies of all particles in the lab-system are greater than 
w 4 A MeV already for the collision at 17 A MeV. On 
the contrary, the solid angle coverage in backward di- 
rection is limited. It has to be noted that phase space 
cells with small average efficiency have large correction 
factors, which amplify the statistical fluctuations of the 
experimental data, and cells with zero efficiency cannot 
be corrected at all. This is the reason why in the back- 
ward hemisphere of the laboratory reference frame (cor- 
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FIG. 2: Average total charge Z to t of forward-going products 
(full dots) and average secondary charge Z scc of PLF (open 
dots) in 93 Nb + 93 Nb and 116 Sn + 116 Sn at 30 A MeV. Dashed 
lines indicate the value of the projectile charge. 

responding in cm. to vu < — 43mm/ns for the reac- 
tion at 38 A MeV and to Uii < — 33mm/ns for that at 
23 A MeV) the applied corrections are not as effective as 
in the forward hemisphere. 

Finally, as the results presented in this paper concern 
LCPs and IMFs emitted in peripheral binary reactions, 
the data need to be corrected for the presence of a back- 
ground of events with a higher multiplicity of heavy frag- 
ments, of which only two have been detected. This back- 
ground of incomplete events (a few percent) has been es- 
timated using the measured three-body events (with the 
procedure described in Ref. |2^) and subtracted from all 
the data presented in this paper. 

D. The particle multiplicities 
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FIG. 3: (color online) Experimental efficiency-corrected av- 
erage multiplicities of (forward-emitted) p, d, t, a and IMFs 
(Z=3-7) as a function of TKEL, for 93 Nb+ 93 Nb (full sym- 
bols) and 116 Sn+ 116 Sn (open symbols). The bombarding en- 
ergies are 17 (triangles, dotted lines), 23 (squares, dot-dashed 
lines), 30 (stars, dashed lines) and 38 A MeV (circles, full 
lines) ; the lines are to guide the eye. 



The average multiplicities of charged particles were 
obtained from the (efficiency corrected and background 
subtracted) experimental distributions of p, d, t, a and 
IMFs, some examples of which are shown in Fig. ^ 

An advantage of using symmetric collisions is that the 
forward-going particles (those with v\\ > in cm.) must 
have the same average characteristics as the backward- 
going ones. Therefore, because of the already mentioned 
much better quality of the data, all presented multiplic- 
ities refer only to particles emitted in the forward hemi- 
sphere of the cm. reference frame (of course, the aver- 
age multiplicities for the whole colliding system can be 
obtained by simply doubling the presented values). 

The symmetry of the system also allows to check the 
quality of the applied efficiency corrections. In fact, 
whatever the reaction mechanism, by adding up the 
charges of forward emitted LCPs and IMFs to the charge 
of the PLF residue, one should reproduce -on average- 



just the projectile charge (of course, this can be done only 
for that subset of the events in which the PLF residue 
was detected and identified in one of the Silicon tele- 
scopes). As an example of the quality of this check, the 
full dots of Fig. |21 show the average total forward-emitted 
charge for the 93 Nb + 93 Nb and 116 Sn + 116 Sn collisions 
at 30 A MeV: within about half a charge unit, the value of 
the projectile (dotted line) is well reproduced, thus giving 
confidence in the present analysis. The rapidly widening 
gap between the full dots and the open ones (which rep- 
resent the average secondary charge Z sec of the detected 
PLF) gives a visual indication of the rising amount of 
charges removed by the LCP and IMF emissions with 
increasing TKEL. 

The obtained average multiplicities (of forward-going 
particles) are shown in Fig. [21 for the systems 93 Nb + 
93 Nb (full symbols) and 116 Sn + 116 Sn (open symbols) 
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at the bombarding energies of 17, 23, 30 and 38AMeV 
(triangles, squares, stars and circles, respectively). The 
statistical errors are smaller than the symbols and the 
lines are just to guide the eye. 

Although data corresponding to more violent collisions 
have been acquired too, the present analysis is limited 
to peripheral and semipcriphcral collisions through a re- 
striction of the TKEL range (from <600 MeV for the 
reaction at 17 A MeV up to <1000 MeV for those at 
38 A MeV). This has the twofold advantage of allowing 
a clear and unambiguous distinction between PLF and 
TLF (in cm. the forward-flying heavy remnant is always 
the PLF) and of limiting the study to regions where the 
binary exit channel is the dominant one. In the con- 
sidered TKEL ranges the binary exit channel remains 
approximately symmetric in mass, as indicated by com- 
parable values of the cm. velocities of the two main frag- 
ments. Only in the last considered TKEL bins weak tails 
of more asymmetric mass splittings appear. In the analy- 
sis, the requirement 0.4 < VplfY^plf + Vtlf) — 0-6 nas 
been used to reject these tails, which are strongly con- 
taminated by incompletely detected three-body events. 

Qualitatively, all multiplicities display a similar behav- 
ior, with a strong dependence on TKEL and a much 
weaker one on bombarding energy and mass of the sys- 
tem. Starting from small TKEL they all increase rapidly 
but then tend to flatten at large TKEL values. Over 
the range of TKEL considered here, the multiplicities 
span about one order of magnitude -or slightly more- for 
the most abundant species (protons and a particles, with 
multiplicities up to several particles per event) and about 
two orders of magnitude for the rarer reaction products 
(which barely reach multiplicities of one per event at the 
highest values of TKEL). 



E. Evaporative and midvelocity emissions 

In order to disentangle the midvelocity emission from 
the sequential evaporation of the PLF, it is convenient 
to make a coordinate transformation into the reference 
frame of the PLF, namely a frame with the i>||-axis still 
oriented along the asymptotic PLF-TLF separation axis, 
but with the origin on the PLF itself. A relativistic 
Lorentz transformation has been applied (instead of sim- 
pler Galilean one), in order to avoid distortions of the 
angular distribution of the fastest particles; this comes 
out to be necessary in particular for protons, which may 
have lab- velocities as large as 30% of the speed of light. 
An example of the obtained yield in the (v PLF , v J LF ) 
plane is shown in Fig. 0fa) for protons, at TKEL= 800 
MeV, in the reaction 93 Nb + 93 Nb at 38 A MeV. The cor- 
responding angular distribution, presented in Fig. 0Jb), 
is obtained by plotting the proton yield as a function 
of cos(6'plf), where #plf is the polar emission angle be- 
tween the velocity of the protons in the PLF frame and 
the PLF-TLF separation axis, as sketched in Fig. 0Je): 
thus cos(#plf) = 1 corresponds to a forward emission 




FIG. 4: (color) Experimental (a) and simulated (c) yields of 
protons in the (ti[' LF ,iJi lF ) plane (with origin in the PLF ref- 
erence frame) and corresponding experimental (b) and simu- 
lated (d) distributions of dcr/dcos(6lpLF) for TKEL=800 MeV 
in the reaction 93 Nb+ 93 Nb at 38 A MeV. Only forward-going 
particles (i.e., with v« in the lab-system larger than the cm. 
velocity) are considered -as sketched in (e)- and this causes 
the yields to vanish for v ii < -35 mm/ns in panels (a) and (c). 



along the flight direction of the PLF in the cm. system. 

For comparison, Fig. Hfc) and (d) show the corre- 
sponding results obtained with a Monte Carlo simula- 
tion, where only an isotropic evaporation from the fully 
accelerated PLF and TLF is considered, as sketched by 
the two circles in Fig. 0Jc). At small TKEL the emis- 
sions from PLF and TLF are well separated from each 
other and the cos (#pLF)-distribution is flat. At larger 
TKEL, on the contrary, the distribution presents a dip 
at cos(#plf)= -1 and a bump at less negative values of 
cos(#plf), as shown in Fig.|3Jd). This distortion of the 
backward part of the angular distribution is due to the 
merging of the two evaporative emissions (from PLF and 
TLF) when the relative velocity of the two sources is re- 
duced: as shown in Fig. E^e), selecting the forward-going 
particles (in cm.) one misses the dash-dotted part of 
the PLF circle, but includes the dashed part of the TLF 
circle. Nonetheless, thanks to the symmetry of the sys- 
tem, the yields and the multiplicities are correctly eval- 
uated. In fact, with respect to the flat behavior of the 
forward part of the cos (£?pLF)-distribution [dashed line 
in Fig. Efd)] , the excess area in the bump neatly com- 
pensates the deficit around -1. 

When comparing the experimental angular distribu- 
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tion of Fig. QJb) with the simulated one of Fig. 0{d), it 
becomes evident that: (i) The shapes are qualitatively 
similar, (ii) The forward part of the experimental angu- 
lar distribution, in the vicinity of cos (#plf) = 1 is indeed 
flat, as expected for an evaporative (nearly isotropic) 
component, (iii) The broad bump around cos(#plf) = 
- 0.5 overcompensates the dip near -1, thus confirming 
the presence of an additional backward emission (with 
respect to the PLF frame), namely an emission in the 
midvelocity region of the phase-space, (iv) The tail of 
the bump extends well into the forward hemisphere in 
the PLF frame 0. 

Thus, as already noted in the analysis of the experi- 
mental angular distribution suggests the superposition of 
two emission components, namely one from the PLF and 
one from "midvelocity" . As it will be shown in Sec. IIV Bl 
there are good arguments for interpreting the first com- 
ponent as an evaporation from the excited PLF; there- 
fore the two terms, evaporation and PLF-emission, will 
be used interchangeably in the rest of this paper. On the 
basis of the previous point (ii), a decomposition can be 
attempted, provided that the shape of the whole evapo- 
rative component can be reliably estimated. Since evapo- 
ration must be forward-backward symmetric in the frame 
of the emitter, it is usual (see, e.g., 0>|2!|) to estimate the 
total evaporation from the PLF by simply doubling its 
forward emission. This may work at higher bombarding 
energies (where evaporation from PLF/TLF and midve- 
locity emissions are well separated in phase space), but 
the above point (iv) shows that, even at the highest en- 
ergy of 38 A McV, the usual procedure would result in a 
significant contamination from the midvelocity emission. 
Which part of the angular distribution can be considered 
sufficiently clean depends on the bombarding energy, on 
the chosen TKEL bin and also on the considered parti- 
cle species. However, for the reactions of this paper, the 
range 0° < 6» PLF < 45° (i.e., 0.7 < cos(6» PLF ) < 1) is a 
reasonable compromise which can be used in almost all 
cases. (We verified that the total evaporative multiplic- 
ities remain stable within ±10% if the data are taken, 
e.g., in the narrower range 0° < 6>plf < 30° while pro- 
gressively larger deviations appear if data beyond 45° are 
included, see Appendix B.). 

In order to estimate, from the data measured at for- 
ward angles, the total yield of the PLF component, it is 
necessary to make some hypothesis on the angular distri- 
bution dcr/d cos(6*plf), which depends on the spin of the 
emitter. In fact this distribution is flat only in case of 
an isotropic evaporation from a zero-spin source, while 
in case of non-zero spin the evaporation tends to con- 
centrate in a plane perpendicular to the spin direction, 
giving origin to an U-shaped cos(#plf) distribution: the 
larger the spin (or the heavier the evaporated particle), 
the stronger the effect. According to a detailed study of 
the correlations of the emitted particles with both the 
PLF and TLF jH^, it was estimated that in the re- 
action 93 Nb + 93 Nb at 38 AMeV the spin of the PLF is 
negligible at low TKEL, but rises to 15±5fi and 30 ± 10ft 



at TKEL w 500 and 800 MeV, respectively. However, the 
shape of the experimental distribution is not very sensi- 
tive to the spin value. In fact, as verified with Monte 
Carlo simulations, even in case of an initial strong align- 
ment of the spin perpendicular to the reaction plane, the 
anisotropy of the experimental angular distribution is re- 
duced - with respect to the theoretical one [^j, |3^ - by 
the misalignment of the spin (caused by the particle evap- 
oration) and by the fluctuations of the reconstructed re- 
action plane. Both effects increase with increasing num- 
ber of particles emitted along the evaporation chain and 
hence with increasing TKEL. Therefore the yields for the 
evaporative component have been obtained by assuming 
a flat distribution (estimated from the data in the for- 
ward range cos(#plf) > 0.7) and applying a correction 
(estimated with Monte Carlo methods) which is at most 
of the order of 15%. 

Once an estimate of the emission of particles from 
the PLF has been obtained, the midvelocity component 
might be derived from the total multiplicities of Fig. fus- 
ing a subtraction procedure Q- However, for an unbiased 
determination of the evaporative component one should 
take into account that, due to momentum conservation, 
the emission of particles - especially of heavy ones - from 
a source of finite mass produces recoil effects which per- 
turb the velocity of the source itself. Altogether, the net 
result of a chain of evaporation steps is a superposition 
of uncorrclatcd perturbations leading just to a smearing 
of the original source velocity, with an average null effect. 
However if the heavy remnant is observed in coincidence 
with a particle emitted in a specific direction, this fact 
introduces a correlation between the velocities of the two 
objects. In particular, if one selects events characterized 
by the emission of certain particles in a restricted angular 
range in the source frame, in those events the perturba- 
tion of the source velocity will have a non-zero average 
value. 

The requirement that a given particle be emitted in 
forward direction selects heavy residues which systemat- 
ically recoil in backward direction. This causes a sys- 
tematic decrease of the cm. velocity of the residue (and 
hence an ovcrcstimation of TKEL), while all the other 
unmeasured or unselected particles contribute only to 
the smearing of the data. As a consequence, the PLF 
multiplicities evaluated in a given bin of TKEL do actu- 
ally pertain to a range of somewhat smaller TKEL values. 
Or, in other words, they erroneously include some contri- 
bution from events belonging to the previous bin (lower 
TKEL) and miss some contribution from events which 
are classified in the following one (larger TKEL). It is 
worth noting that the heavy residue recoils in the (op- 
posite) forward direction when it evaporates particles at 
angles 6*plf > 90°, which then merge with the midveloc- 
ity emissions. In this case the cm. velocity of the residue 
increases systematically and the resulting TKEL is un- 
derestimated. This has to be taken properly into account 
before subtracting the evaporative component from the 
total multiplicities. 
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FIG. 5: (color online) Experimental, efficiency-corrected evaporative (left) and midvelocity (right) multiplicities of forward- 
emitted p, d, t, a and IMFs (7=3-7) as a function of TKEL, for 93 Nb+ 93 Nb (full symbols) and 116 Sn+ 116 Sn (open symbols). 
Error-bars indicate not statistical errors (they are smaller than the symbols), but rather variations of multiplicity values with 
analysis procedure (see text). The symbols, with the same meaning as in Fig. are slightly displaced horizontally in order to 
appreciate the individual error bars at a given TKEL. 



IV. RESULTS 
A. Particle multiplicities 

The evaporative (Alovap) and midvelocity (.M m idv) 
multiplicities of forward-going particles are presented, 
as a function of TKEL, in the left and right part of 
Fig. UJ respectively. The multiplicities refer -with the 
same symbols- to the same experimental data used for 
Fig. |21 The vertical scales differ for different particles, 
but they are the same for the evaporative and midveloc- 
ity components of a given particle. It is also worth noting 
that at 17 A MeV it is not possible to reliably extract the 
midvelocity multiplicities Al m idv for protons and a par- 
ticles, because of the overwhelming contribution of the 
evaporative part. 

The error bars of Fig. UJ do not represent statistical 
errors (which are usually smaller than the symbols). In- 



stead, at each point, they show the maximum variation 
obtained by choosing a different angular window to select 
the evaporative component, or by changing the extrap- 
olation to obtain the whole angular distribution, or by 
switching on/off the recoil corrections. More details are 
given in Appendix B. 

Let us first consider the average multiplicities for the 
evaporative emission, A^ C vap- They all display a mono- 
tonic increase with increasing TKEL: steeper for the most 
peripheral collisions, flatter for the less peripheral ones. 
However, the magnitude of this increase is different for 
the different particle species. In fact, at small TKEL, the 
average multiplicities range from a few tenths for protons 
down to about 10~ 3 for tritons or even 10 -4 for IMFs, 
while at the highest TKEL values considered here the 
multiplicities are much more leveled off, with 2-3 pro- 
tons or a particles, as compared to 0.2-0.3 tritons or 
IMFs. Thus, the entire evolution with TKEL may be 
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comprised within about a decade (protons in the Nb+Nb 
system), or it may span more than three orders of mag- 
nitude (IMFs). 

Concerning the dependence on the system size, Fig. [3] 
shows that, at a given TKEL, the evaporative multiplici- 
ties in the 116 Sn+ 116 Sn collision (open symbols) are sys- 
tematically s mailer than the corresponding multiplicities 
in the 93 Nb+ 93 Nb collision (full symbols). This holds 
true for all particle species and at all TKEL values, but 
it is more clearly visible in the panels for a particles and 
protons and generally in the lowest TKEL bins. This ef- 
fect is well reproduced by calculations with the statistical 
code Gemini and may be explained by the fact that 
nuclei produced in the 116 Sn + 116 Sn reaction, because of 
the slightly larger initial value of N/Z (1.32 versus 1.27 
of 93Nb), get rid of their excitation energy by emitting 
more neutrons and less charged particles. 

As for the beam-energy dependence, at a given TKEL 
there is a limited sensitivity of A^evap to this parameter, 
but nonetheless one can clearly observe, for all particle 
species, the tendency to rise with decreasing beam en- 
ergy, a feature that becomes more prominent for heavier 
particles. For example, in Nb+Nb, the multiplicities are 
nearly the same for protons at all four bombarding ener- 
gies, while for a particles and IMFs at small TKEL they 
rise by about a factor of 1.5-2 from one beam-energy to 
the next lower one. A possible explanation is that, as 
already noted, with increasing beam energy the midve- 
locity emissions draw off an increasing amount of mass 
and energy from the collision, so that a given TKEL cor- 
responds to smaller excitations of the primary PLFs and 
hence to lower values of A4 evap . Moreover, at high bom- 
barding energies the contact time is shorter and the an- 
gular momentum transfer is likely to be smaller so that 
- in spite of the larger values of the orbital angular mo- 
menta in the entrance channel - the primary PLFs may 
reach lower spin values, an additional fact hindering the 
evaporation of heavier particles. 

In any case, the variations of M cvap with bombard- 
ing energy are small when compared with the variations 
(about two orders of magnitude, or more) occurring over 
the full range of TKEL. Thus, at Fermi energies the vari- 
able "TKEL" is still reasonably well correlated with the 
total amount of particles emitted from the PLF. This 
observation supports the conclusion j23| that, at these 
beam energies, equal values of TKEL still indicate com- 
parable excitation energies deposited in the PLF (and in 
the TLF as well). 

The midvelocity multiplicities Al m idv, displayed in the 
right part of Fig. present an increase with TKEL simi- 
lar, at first sight, to that of the evaporative multiplicities 
■Afevap- However, a closer inspection shows that their 
behavior is in many respects different and to a certain 
extent complementary. First of all there is a less pro- 
nounced rise, the entire evolution with TKEL spanning 
-for all particles- less than two decades; in some cases the 
multiplicities seem even to decrease at the highest TKEL. 
Then, the dependence on the mass size favors the heav- 
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FIG. 6: (color online) Ratios of midvelocity-to-total multi- 
plicities of forward-emitted p, d, t, a and IMFs, as a function 
of TKEL. Symbols have the same meaning as in Fig. [3] 



ier 116 Sn+ 116 Sn system, where, at a given TKEL, one 
generally observes larger A^midv values, a feature par- 
ticularly evident for hydrogen isotopes. Finally, at vari- 
ance with what is observed for .Mcvap, at a given TKEL 
value the midvelocity multiplicities show an appreciable 
increase with increasing beam energy, especially for what 
concerns the lighter particles. 

In order to better appreciate the different dependence 
of the evaporative and midvelocity multiplicities on the 
beam energy and TKEL, Fig. HJ presents the ratios of 
the midvelocity components to the total multiplicities 
measured in the collisions 93 Nb + 93 Nb and 116 Sn + 
116 Sn. One observes that generally, for a given beam en- 
ergy, all ratios decrease with increasing TKEL. Moreover, 
with decreasing beam energy, the ratios have smaller val- 
ues and display a faster decrease with increasing TKEL. 
With the notable exception of the protons and - to a 
lesser extent - of the a particles, at low TKEL (i.e., in 
the most peripheral events) the multiplicity ratios tend 
to reach values as large as 0.8-0.9. This fact confirms 
that peripheral collisions are the best environment to in- 
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FIG. 7: Dependence on beam energy of evaporative (up- 
per panels) and midvelocity multiplicities (lower panels) 
for forward-emitted p, d, t, a and IMFs (Z=3-7) in the 
93 Nb+ 93 Nb collision. The data correspond to TKEL values 
of 200 and 400 MeV (left and right panels, respectively). 
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FIG. 8: (color online) Multiplicities of PLF-emitted parti- 
cles as a function of l/^/TKEL/2. for 93 Nb+ 93 Nb at 17, 
23, 30 and 38 A MeV. Lines are linear fits to the logarithm 
of the multiplicities in the low-energy part of the excitation 
functions. 



vcstigate the phenomenon of the midvelocity emissions 
with the least contamination from sequential evapora- 
tion. This prevalence of the midvelocity particles de- 
creases with increasing centrality of the collision. Com- 
paring the different panels, one sees that protons are 
mainly produced in evaporative processes and a particles 
are almost evenly shared between the two mechanisms. 
On the contrary, the most "exotic" products, i.e., tritons 
and especially the IMFs, are those which are most specific 
of the midvelocity emissions: at 38 A MeV, even at the 
highest TKEL, more than 70% of all emitted IMFs are 
attributed to the midvelocity component and this per- 
centage rises to almost 100% at small TKEL. 

The different evolution of M evap and M m idv with bom- 
barding energy is shown in Fig. 0for two TKEL bins in 
the collision 93 Nb + 93 Nb. While all evaporative multi- 
plicities stay almost constant, or even show a weak de- 
crease with increasing bombarding energy, the midveloc- 
ity multiplicities display a general increasing behavior, 
more pronounced for the data at higher TKEL. 



B. Nature of the emissions 

In spite of repeated experimental and theoretical ef- 
forts over the years, the question about the nature of the 
mechanism(s) responsible for the observed LCP and IMF 
emission is still not completely settled. 

For what concerns the PLF emissions, it is natural 
to expect that an evaporation-like de-excitation should 
be present. Indeed the usual procedure adopted in the 
literature (see, e.g., Q), and also in the present paper, 



for separating the two emission components is based on 
that working hypothesis. Its validity can be checked by 
looking at several characteristics of this emission, in com- 
parison with the results of calculations with a statistical 
evaporation code like Gemini [33| ■ Hereafter we present 
some pieces of information which support the hypothesis. 

According to statistical theory, the partial decay width 
associated with a given exit channel should present an ex- 
ponential dependence on the temperature T of the source 
(Boltzmann factor, mainly due to the increase of level 
density with increasing T) of the type: T oc exp(— B/T), 
where B is the barrier associated with that channel. 
Therefore also the average multiplicities M. of the various 
evaporated particles may be expected to present a similar 
exponential dependence on 1 /T, or - in a Fermi gas model 
- on the inverse square root of the excitation energy, 
M. cx cxp(— c/vE*), where c is a constant, dependent 
on the particle species. In order to compare data with 
these expectations, the first non-trivial task is to estimate 
the PLF excitation energy. At low bombarding energies 
(where midvelocity emissions are negligible and TKEL 
represents the total excitation of the system) and for sym- 
metric systems, TKEL/2 is certainly a good average esti- 
mate of the excitation energy of each of the two primary 
products of a binary reaction. As already noted, with in- 
creasing bombarding energy this interpretation becomes 
more and more questionable due to the increasing rele- 
vance of the midvelocity emissions. However, the overall 
energy balance of the reaction 93 Nb+ 93 Nb at 38 A MeV 
performed in Ref. [f| has demonstrated the persistence of 
an approximately linear correlation between TKEL and 
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FIG. 9: (color online) Average excitation energy of PLF 
(left) and midvelocity "source" (right) as a function of TKEL 
for 93 Nb+ 93 Nb at 38 A MeV. Uncertainties due to different 
analyses and assumptions about the n-richness of the midve- 
locity emissions are indicated by error bars. Dashed lines and 
shaded areas show estimates and uncertainties of Ref. 0. 
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FIG. 10: (color online) (a) Multiplicities of PLF-emitted par- 
ticles as a function of l/y^E*^, for 93 Nb+ 93 Nb at 38 A MeV, 
with Ep LF excitation energy (in MeV) of the primary PLF 
(symbols as in Fig. [HJ; lines are linear fits to the points, (b) 
Similar plot, for the decay of a 93 Nb nucleus according to 
Gemini calculations. 



the excitation energy of the two primary products. Thus 
one can suppose that TKEL/2 still represents a measure 
of the excitation energy of the PLF (and of the TLF as 
well), although it may be that the scale is no longer a 
quantitative one, differing at the various bombarding en- 
ergies. In Fig. |SJ the data of the reaction 93 Nb+ 93 Nb at 
the four bombarding energies of 17, 23, 30 and 38 A MeV 
are presented |22l | in a semilogarithmic plot as a function 
of 1/ ^/TKEL/2. Indeed, the logarithms of the multi- 
plicities of all particles display the same behavior at all 
bombarding energies, namely a linear one. This fact is 
put into evidence by the lines, which arc fits to the points 
in the low TKEL region, for each particle species. 

A more quantitative scale of excitation energy can be 
estimated following the analysis performed in Ref. Q 
for the reaction 93 Nb+ 93 Nb at 38 A MeV. There the av- 
erage excitation energy of the primary PLF and of a hy- 
pothetical midvelocity "source" was obtained from the 
total energy balance of the reaction, i.e., by summing 
up the measured kinetic energies of the respective emit- 
ted particles and taking into account the average Q- value 
for disassembling the system into the final reaction prod- 
ucts. For the semipcriphcral events of 93 Nb+ 93 Nb at 
38 A MeV, the excitation energies of the two sources have 
been re-evaluated with some improvements in the anal- 
ysis and the results are shown in Fig. The previous 
estimates of @ (dashed lines) are indicated too, together 
with their uncertainties (shaded areas) resulting from two 
rather extreme assumptions about the n-richness of the 
midvelocity emissions (from N/Z= 1.1 to 1.44). For the 
present evaluation, the uncertainties due to the same hy- 
potheses on N/Z and to possible variations in the analysis 
(see Appendix B) are indicated by error bars. Most of 
the difference with respect to the previous results of Q 
is due to the use of relativistic kinematics and to the 
correction for recoil effects. 

Figure iFOT a) shows the PLF particle multiplicities as 
a function of the inverse square root of the excitation 
energy of the primary PLF, 1 / y^Ep LF . [The first point 
(TKEL = 100 MeV) is not used here because of the large 



relative uncertainty of its position on the x-axis, depend- 
ing on the analysis method (see Appendix B).] When 
this more appropriate estimate of the excitation energy 
of the "source" is used, an improved linear behavior is 
apparent in the logarithmic presentation, as shown by 
the displayed linear fits to the points. Thus the observed 
shape of the excitation functions is compatible with a 
process ruled by statistical laws. 

A similar plot, based on the multiplicities calculated 
with the code Gemini for the decay of an excited 93 Nb 
of spin 30 H, is shown in Fig. HOf b). The agreement of 
the calculations with the experimental data for the PLF- 
emissions is rather good, thus confirming the evaporative 
origin of this component. No attempt has been made to 
take into account, in a detailed way, possible changes of 
mass and spin of the primary PLF with increasing E*, 
as this was beyond the scope of the comparison. How- 
ever the relatively weak sensitivity of the calculations to 
these parameters and the good agreement with the ex- 
perimental data for the PLF emissions are demonstrated 
in Fig. 1111 In this picture the open symbols show the 
multiplicities of n, p, d, t, a particles and light IMFs 
(Z=3-7) calculated with the Gemini code for two val- 
ues (E*= 100 and 200 MeV) of the excitation energy of 
the evaporating source and for different assumptions on 
its mass and spin. The full dots show the corresponding 
experimental multiplicities of PLF-emitted charged par- 
ticles (neutrons are not measured) for two TKEL bins 
which, according to the energy balance of Fig. EJa), cor- 
respond to the selected values of E*. Again, the good 
agreement between experimental data and calculations 
indicates that the decay of PLFs produced in peripheral 
and semipcriphcral collisions is compatible with the usual 
evaporation of an excited nucleus at normal density, at 
least for what concerns the first 3-4 fm of overlap (see 
Appendix A). So, at variance with the claims of other au- 
thors Q , in order to reproduce our data there is no need 
to resort to models (like SMM [34[) which describe the 
fragmentation of an expanding diluted nucleus (typically 
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FIG. 11: (color online) Multiplicities of light particles and 
light IMFs (Z=3-7) evaporated from the PLF in 93 Nb+ 93 Nb 
at 38AMeV, for an estimated PLF excitation energy of 100 
and 200 MeV (left and right panel, respectively) . Experimen- 
tal data (full dots) are compared with Gemini calculations 
(open symbols) with different mass and spin of the source; 
for the calculations also the neutron multiplicities are shown. 



at one third of normal density). Possibly, the reason is 
that the events analyzed by the cited authors are less pe- 
ripheral (indeed their estimated PLF excitation energies 
are around 4-5 A MeV, i.e., higher than ours) and cor- 
respond to a rather particular event selection (at least 
3 IMFs) and not to the bulk of the collisions as in our 
case. It has also been verified that, as expected for an 
equilibrated emission, the velocity spectra of the emitted 
LCPs have approximately Maxwellian shapes with slope 
parameters of the order of 2-5 MeV, depending on the 
selected TKEL value. 

Finally, the lower part of Fig. ^| shows the average 
N/Z ratio of hydrogen isotopes emitted by the PLF in 
the reaction 93 Nb+ 93 Nb at 38 A MeV (full dots). This 
ratio increases with TKEL from 0.1 to about 0.4 and is 
compatible with the results of Gemini calculations for 
the decay of a 93 Nb nucleus at normal density with the 
appropriate excitation energy obtained from Fig. [ffia). 

For what concerns the midvelocity component, some 
characteristic aspects (like, e.g., the emission pattern of 
LCPs and IMFs, the space-time extension of the "source" 
and its associated energy 

U E II 13 El ) have been al- 
ready investigated in the past. However, very different 
hypotheses have been proposed for the underlying mech- 
anism, ranging from fully dynamical processes (like, e.g., 
surface instabilities in the non-spherical, transient shapes 
of the interacting system to purely statistical ones 
(like, e.g., a proximity-enhanced statistical decay of PLF 
and TLF in the external inhomogeneous Coulomb field 
of the other reaction partner although it cannot 

be excluded that more than one mechanism contributes 
to the observed phenomena Q . 

The average N/Z ratio of the hydrogen isotopes in the 
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FIG. 12: Average value of N/Z for evaporative (full dots) and 
midvelocity (open dots) emission of Hydrogen as a function 
of TKEL in 93 Nb+ 93 Nb at 38 A MeV. Dashed lines are the 
results of calculations with the code Gemini. 



midvelocity component is shown by the open dots in the 
upper part of Fig. ^1 again for the reaction 93 Nb+ 93 Nb 
at 38 A MeV. This ratio has now much larger values (from 
0.6 to about 1.0), thus indicating a substantial difference 
between the emissions from the PLF and those at midve- 
locity. This may be an indication of neutron enrichment 
at midvelocity, as proposed on the basis of the emission 
of complex particles Q, 0, 0, LH Ell IHh or it mav be 
somehow related to the reduced size of the "source" and 
hence to its higher energy concentration, as indicated in 
□ 

The energy deposited in the matter at midvelocity, 
-^midv estimated from the energy balance of Fig. OJb), is 
used in Fig. 1131 to present also the midvelocity multiplic- 
ities 2 in a logarithmic plot as a function of 1 / y/E^ idv : 
to our knowledge, this is the first time that midveloc- 
ity multiplicities are presented in this way. (Here too, 
the first point at TKEL = 100 MeV has not been used.) 
First of all, we want to point out some remarkable dif- 
ferences between the midvelocity emissions of Fig. 1131 
and the PLF emissions of Fig. HOf a). For example, all 
midvelocity multiplicities are compressed within about 
one decade (note the different horizontal and especially 
vertical scales of the two figures) and consequently their 
slopes are all sizably flatter. Even more striking is the 
fact that the relative abundances of the emitted particles 
are different, with remarkable inversions between protons 
and a particles and between tritons and IMFs. 

In this presentation as a function of 1/ \/E*^ idv , one ob- 
serves a linear correlation also for the midvelocity emis- 
sions, as demonstrated by the linear fits to the points, 



2 Here, as in Ref. the particle multiplicities and ^^idv of the 
midvelocity "source" refer to the forward-going particles only; 
for the total midvelocity processes one has to double both the 
multiplicities and the estimated E* . , . 

1 miav 
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FIG. 13: (color online) Semilogarithmic plot of midvelocity 
particle multiplicities as a function of l/v^S^ idv for Nb+Nb 
at 38AMeV (symbols as in Fig. E^ idv is the excitation 
energy of the emitting "source", lines are fits to the points. 



but no easy interpretation is at hand. For a different 
process, namely multifragmentation in central collisions, 
it was argued [3i| that such a linear behavior (called 
thermal scaling) would be the proof of the statistical, se- 
quential decay of the "source" (the role of dynamics being 
relegated to its formation), but this conclusion is amply 
disputed. In fact, very different models proved able to 
reproduce a similar behavior, and no general consensus 
is reached on aspects like the appropriate treatment of 
the data and estimation of the source temperature, the 
sequentiality or simultaneity of the process , th e role of 
fluctuations and correlations (see, e.g., [321 iM LlH and 
references therein). Although the linear behavior in itself 
may be suggestive of a thermal process, its mere obser- 
vation is not a proof and one actually expects that in 
peripheral collisions the dynamics should play a much 
more important role than in central collisions. 

As a last point, it is worth noting that these striking 
differences between evaporative and midvelocity multi- 
plicities give strong support to and demonstrate the ef- 
fectiveness of the separation procedure outlined in Sec. 
IIII El which has been used to disentangle the two com- 
ponents. 



V. SUMMARY AND CONCLUSIONS 

The symmetric collisions 93 Nb + 93 Nb at 17, 23, 30 
and 38AMeV and 116 Sn+ 116 Sn at 30 and 38AMeV 
have been investigated using the FIASCO setup, a low- 
threshold multidetector covering a large fraction of the 
forward solid angle. 

The present analysis is focused on peripheral and 
scmipcripheral events resulting from binary or quasi- 
binary reaction processes. They are characterized by the 
presence in the exit channel of two heavy remnants, which 
can be identified as the projectile- and target-like frag- 
ments. These binary or quasi-binary processes are found 
to exhaust more than half of the expected total reaction 
cross section at all the investigated beam energies. 



The kinematic variable TKEL has been used for select- 
ing event samples corresponding, on average, to decreas- 
ing values of the impact parameter. It is worth stressing 
that TKEL, which at low bombarding energies represents 
a good estimate of the total excitation of the colliding nu- 
clei, looses such a physical meaning at Fermi energies and 
therefore it is used here only as an ordering parameter. 

The average multiplicities Ai of light charged parti- 
cles and intermediate mass fragments with 3 < Z < 7 
have been obtained as a function of TKEL for all the in- 
vestigated systems. The data have been evaluated in a 
range of TKEL values which is estimated to roughly cor- 
respond to the outermost 30% of the impact parameter 
range leading to nuclear interaction. The emission pat- 
tern of these particles in the (v ii ,i>j_) plane - the velocity 
components being referred to the PLF - TLF separation 
axis in the cm. reference system - clearly shows the 
presence of two components, one representing an emis- 
sion from the excited PLF (and TLF) and a second one 
at midvelocity, i.e., at velocities intermediate between 
that of PLF and TLF. A careful analysis of the shape 
of the angular distribution in the reference frame of the 
PLF has been applied in order to distinguish and esti- 
mate the multiplicities of particles emitted by the PLF 
and, by a subtraction procedure, those of the particles 
emitted in the midvelocity region. Correction have been 
applied to take into account the efficiency of the setup, 
as well as physical effects due to spin and recoil of the 
emitting nuclei. 

Both emission components increase with increasing 
TKEL (i.e., decreasing impact parameter), although 
their behavior is different. More "exotic" particles (like 
tritons, IMFs and to a lesser extent deuterons) are char- 
acteristic of midvelocity processes in the most peripheral 
events, where they outnumber the emissions from PLF. 
Moreover, for a fixed TKEL, the midvelocity emissions 
present a clear increase with increasing beam energy, 
while the emissions from PLF show little dependence on 
the beam energy. 

The dependence of PLF emissions on the estimated 
excitation energy of the PLF follows that expected for a 
decay mechanism governed by a barrier and is well de- 
scribed within statistical models. In fact, calculations 
with the statistical code Gemini well reproduce all ex- 
perimental features, including the slopes of the excitation 
function, the relative and absolute abundances of the var- 
ious particle species and also the isotopic composition 
of Hydrogen particles. Thus, at least for the peripheral 
and semiperipheral collisions, the decay of the PLF (and 
TLF) is in many aspects compatible with the usual evap- 
oration from an excited nucleus at normal density. 

Surprisingly, also the midvelocity emissions display 
a similar type of dependence on the amount of en- 
ergy which is localized in the hypothetical midvelocity 
"source". However, this mechanism is quite different 
from the usual evaporation, as it is shown by the inver- 
sion in the relative abundance of the emitted particles 
and by the tendency to emit more n-rich light particles. 
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Indeed, the emission pattern observed in this work and 
in previous ones Q, clearly demonstrates that this hy- 
pothetical "source" is not a simple spherical one, sitting 
at midvelocity and isotropically emitting particles; on 
the contrary, in order to reproduce the observed emission 
patterns, one needs a more complex "source" , extended 
both in space and time. For a satisfactory explanation of 
its characteristic features, o ne p robably needs a dynamic 
description of the collision [llj , but its successive decay 
seems to possess some features reminiscent of a statistical 
process. 
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APPENDIX A: IMPACT PARAMETER 
ESTIMATION 

Various experimental variables are used in literature 
for estimating how central or peripheral a collision is 
and for sorting the measured events in (nearly) homo- 
geneous samples. Among the variables related only to 
the main reaction partners, an often used one is the sec- 
ondary charge Z sec of the PLF residue (at low bombard- 
ing energies or for quasi-peripheral events) or the charge 
of the heaviest fragment Z max (at higher energies or for 
more central collisions): in fact, on average, the lighter 
this charge, the more violent and central the collision is 
likely to be |40j |. Such non- kinematic variables may be 
used for studying kinematic aspects of particles (like, e.g., 
their angular distributions), but they are not equally well 
suited for studying multiplicities. Because of the finite 
(and not too large) total charge of the system, spurious 
correlations may appear. 

In order to avoid this problem, it is preferable to sort 
the data according to a kincmatical variable, such as the 
secondary velocity v p^p °f the PLF in the laboratory ref- 
erence system |4l| , or the relative velocity v re i of the two 
main reaction products, or the energy loss per nucleon e* 
of a binary collision [42] ■ In this paper we have adopted 
TKEL, as defined in Eq. QJ. In case of a binary pro- 
cess with frozen initial mass asymmetry one expects to 
find a good correspondence between all these kinematic 
variables, but if the mass asymmetry in the exit channel 
is allowed to change, then TKEL has the advantage of 
explicitly taking into account these variations. 
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FIG. 14: (color) Experimental correlations between the kine- 
matic variable TKEL [see Eq. JTJ] and (a) relative PLF-TLF 
velocity; (b) lab-velocity of the PLF; (c) secondary charge of 
the PLF. Data refer to 93 Nb + 93 Nb at 38 A MeV. 



The experimental correlation of TKEL with v re i and 
v ^ for the data of 93 Nb + 93 Nb at 38 A MeV is shown 
in Fig. I14f a) and (b), respectively. The first correlation, 
TKEL-v re i , is extremely narrow in the whole considered 
range: this is little surprise, since in symmetric collisions 
most binary exit channels have nearly the same reduced 
mass, so that TKEL is approximately proportional to 
■y 2 rel [see Eq. QJ]. The next correlation, TKEL VS V , 
is narrow at low TKEL, but tends to become wider with 
increasing inelasticity of the collision. Finally, Fig. 1141 c) 
also shows the correlation of TKEL with Z scc : the two 
variables are well correlated but with rather large fluctu- 
ations. 

The usage of TKEL as an impact parameter estima- 
tor is supported by the kinematic analysis of events gen- 
erated with the Quantum Molecular Dynamics (QMD) 
code Chimera [43 ■ The velocity vectors of the two main 
fragments produced by such QMD calculations were used 
to deduce the variable TKEL, just with the same kine- 
matic procedure used for the experimental data: the re- 
sult of this analysis shows a nice correlation between im- 
pact parameter and TKEL [2l|. The average value of 
the obtained TKEL, plotted as a function of the impact 
parameter used as input to the calculations, is shown in 
the upper part of Fig. ED for the investigated reactions 
(except for the 17 A MeV case, because at such a low 
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FIG. 15: (color online) Upper panel: average correlation be- 
tween impact parameter and the kinematical variable TKEL 
for events generated by the QMD code Chimera |4^| . Lower 
panel: average correlation between impact parameter and 
TKEL obtained from the integration of the experimental reac- 
tion cross section for 93 Nb+ 93 Nb at 17, 23, 30, 38 A MeV and 
116 Sn+ 116 Sn at 30, 38 A MeV; arrows indicate the cm. avail- 
able energies which are within the displayed energy range. 
Lines are just an aid to guide the eye. 



bombarding energy the QMD approach seems hardly ap- 
plicable). As expected, there is a monotonic increase of 
TKEL when the impact parameter is decreased starting 
from the grazing value 6 graz . Less obvious is the finding 
that, for a given system, at large impact parameters the 
correlations are almost independent of the beam energy. 
The model suggests that, at different bombarding ener- 
gies, peripheral events sampled in the same TKEL bins 
may correspond to collisions with roughly the same im- 
pact parameter. As the collisions become less and less 
peripheral, the curves progressively separate and bend 
down, as they approach - at different TKEL values - 
the limit of their respective available energy in the cm. 
system. Concerning the comparison of the two systems, 
for 116 Sn + 116 Sn the correlations begin at a larger value 
of b, as might have been expected, and proceed almost 
parallel to those of the lighter system 93 Nb + 93 Nb. 

As an alternative, one can perform an association 
of TKEL with impact parameter using a method al- 
ready applied in Ref. [g and based on a refinement 
of an even older recipe (44J. The basic hypothesis is 
that there is a good average correlation between de- 



creasing impact parameter and macroscopic behavior of 
the reaction products in terms of scattering angle and 
kinetic energy dissipation. So one can try to follow 
the evolution of the reaction from the two-dimensional 
plot d 2 a/d(9 c . m .)d(TKEL), the so-called Wilczynski-plot 
(45j , and empirically determine - on average - an impact 
parameter scale by means of an integration of the reac- 
tion cross section, starting from the elastic region, going 
over to the quasi-elastic region and then down into the 
deeply inelastic region; events leading to fusion or to sev- 
eral heavy fragments are assumed to be located at the 
lower end of the impact parameter scale. 

In practice, during our experiment, events from a min- 
imum bias trigger ("singles" , requiring a number M>1 of 
hits in the gas detectors) were acquired at a reduced rate, 
together with rarer events from more selective triggers 
(requiring higher-fold hits). For the four most forward 
gas detectors, "clean" angular distributions can be ob- 
tained for these minimum-bias events by selecting appro- 
priate equal windows in the azimuthal angle and further 
requiring that the time-of-flight be compatible with that 
of elastically scattered projectiles. The so obtained an- 
gular distributions nicely reproduce the l/sin 4 (0 c . m ./2) 
shape expected for the Rutherford scattering of point- 
like charges, until the region near the grazing angle is 
reached. The method is very sensitive; small differences 
in the rates of the four most forward gas detectors, have 
been attributed to a misalignment (generally less than 
one tenth of a degree) of the beam with respect to the 
optical axis and were used to correct the polar angles. 
For each of the investigated reactions, a simultaneous fit 
of the Rutherford formula to the angular distributions 
of the four gas detectors has been used to estimate the 
conversion factor from counts to millibarns, with an un- 
certainty below a few percent. 

The procedure is then illustrated with the help of 
Fig. El where the correlation d 2 a / d9 c . m .<iTKEL for bi- 
nary events is sketched together with its total projection 
der/<iTKEL (on the left) and a projection dcr/d6 c . m . of 
the elastic and quasi-elastic region (on the top). The 
region A is dominated by elastic scattering (which is re- 
sponsible for the large peak around TKEL«0 in the pro- 
jection on the left), while regions B (still at TKEL«0, 
but beyond the grazing angle 6*| r ^ ) and C are popu- 
lated by reactions. The two-body cross sections of Ta- 
ble [H] have been obtained by integrating the efficiency 
corrected experimental yields in regions B and C (with 
a small correction to take into account the quasi-elastic 
cross section lying below the elastic ridge in region A), 
while the three-body cross sections come from the lower 
part of region C only. 

The correlation between impact parameter b and 
TKEL is obtained by filling the triangle-shaped distribu- 
tion dajdb (upper part of Fig. ITfj|> in successive steps. 
Starting from a value of 6 c . m . still in the elastic re- 
gion (which is uniquely related to impact parameter by 
Rutherford scattering), one integrates the measured cross 
section at larger angles (remaining part of the region A 
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FIG. 16: Sketch of the cross section integration 



and region B), thus determining a new impact parameter 
value, which is the first point in the experimental correla- 
tions [see Fig. I15f b1]. Now one proceeds towards smaller 
and smaller impact parameters by integrating the mea- 
sured cross section for more and more inelastic events (in 
region C). In a "sharp cut-off" picture, each successive 
integrated piece of cross section is used to completely fill 
the next (to the left) slice of the triangle da/db, thus 
determining the next point of correspondence between 
impact parameter b and TKEL. Of course, for this pro- 
cedure to work reasonably, what must be integrated is the 
total reaction cross section. For the present work, focused 
on peripheral collisions, this means that it is necessary 
to obtain: 

i) a good quantitative measurement of the most rep- 
resentative exit channels of the reaction (i.e., two- 
body and, to a lesser extent, three-body events); 

ii) a good correction for the inefficiencies of the setup 
( "experimental filter" ) ; 

iii) a smooth transition, without detection gaps, from 
the elastic to the inelastic events (i.e., not too 
complicated triggers and good efficiency for quasi- 
elastic events). 

The results for the investigated reactions are presented 
in the lower panel of Fig. Qualitatively there is 

a good agreement with the results of the QMD esti- 
mates, as the same main features are well reproduced: (1) 
There is the same monotonic decrease of b with increasing 
TKEL. (2) For very peripheral collisions the correlations 
are almost independent of the beam energy, while they 
tend to separate at smaller impact parameters. (3) For 
116 Sn + 116 Sn the curves are displaced by an almost con- 
stant value of about 1 fm to larger impact parameters. 

For each curve, the first point - that at the highest b 
- is the most critical one and it is estimated to be af- 
fected by an overall uncertainty of the order of 1 fm (due 



to corrections of misalignment of the beam axis, fitting 
with Rutherford scattering, integration of elastic cross 
section, matching of elastic and inelastic events). Being 
the curves obtained by successive integrations of the cross 
section, the errors are strongly correlated and the uncer- 
tainty on the first point propagates to all other points; 
in other words, if the first point needed to be moved 
up or down, all the others would move almost rigidly in 
the same direction. Considering that the FIASCO setup 
has been developed and optimized for peripheral reac- 
tions, this impact parameter determination becomes less 
reliable when the collisions become more central In fact, 
one may start missing some relevant exit channel and, in 
turn, this progressive underestimation of the total reac- 
tion cross section may cause the impact parameter to de- 
crease too slowly with TKEL. Indeed, this might explain 
why, at large TKEL and especially for larger beam en- 
ergies, the experimental curves shown in the lower panel 
of Fig. ^)]do n °t bend downwards as much as the QMD 
calculations. 

In the TKEL ranges used for the evaluation of the 
multiplicities presented in this paper (TKEL < 600 and 
800 McV for Nb+Nb at 17 and 23AMeV, respectively; 
<1000 MeV in all remaining cases) the binary channel 
accounts for about 2.0-2.2 barn in the 93 Nb + 93 Nb col- 
lisions and about 1.9 barn in the 116 Sn + 116 Sn collisions 
(corresponding to ~ 40-50% and ~ 35% of (T^, respec- 
tively), with an uncertainty of about 200-300 mb. The 
applied restriction on the ratio of the two cm. velocities 

(0.4 < Uplf/(>pLf + v tlf) < °- 6 ) cuts lcss than 5% of 
the considered cross section and is appreciable only in 
the last considered TKEL bin. From the integration of 
the experimental cross section it is estimated that the 
presented multiplicities refer approximately to the outer- 
most third of the impact parameter range (6 oxp / bgf^ z > 
60-70%). In this region, choosing the same TKEL bin 
allows to compare results at roughly the same impact 
parameter if the system is fixed and only the beam en- 
ergy is varied, while it is confirmed that there is a shift 
of about 1 fm between the estimated impact parameters 
for the two systems 93 Nb + 93 Nb and 116 Sn + 116 Sn. 

APPENDIX B: DEPENDENCE ON THE 
ANALYSIS METHOD 

In order to estimate the uncertainties which may affect 
the determination of the average evaporative and mid- 
velocity multiplicities (-M OV ap and A^ m idv, respectively), 
the data have been analyzed by using slightly differ- 
ent procedures for separating the total multiplicities into 
these two components. More specifically, somewhat dif- 
ferent choices have been taken at three important points 
of the analysis and the induced variations have been con- 
sidered as representative of the sensitivity of the results 
to different evaluation procedures. The three points are 

a) the angular range (in the PLF frame) used to esti- 
mate the evaporative component, 
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b) the successive extrapolation of the angular distri- 
bution to the whole range 0°-180°, 

c) the corrections for recoil effects. 

As for point (a), besides the adopted range of 0°-45°, 
a more conservative one (0°-30°) has been considered, as 
well as the range 0°-90°, which is often used in the litera- 
ture. As already discussed with regard to Fig. this last 
choice pays the advantage of estimating the whole angu- 
lar distribution by means of a simple reflection around 
90° [without the extrapolation of point (b)] with the dis- 
advantage that the evaporative component may be con- 
taminated by the tail of the midvclocity emission. 

The extrapolation at point (b) requires some hypothe- 
sis on the spin value of the evaporating source and its mis- 
alignment during the decay, while for a zero-spin source 
one expects a roughly cos 9 pLF-shaped angular distri- 
bution da/dd plf (although recoil effects may somewhat 
distort it). So the considered alternatives are an analyt- 
ical extrapolation with the cos 9 plf shape and a numer- 
ical one, obtained with a Monte Carlo simulation which 
assumes a non-zero spin value (rising with TKEL from 
h to 40 K) and takes into account the progressive spin 
misalignment along the evaporation chain. 

Figure El sriows the differences which are typically ob- 
tained, by taking as an example the evaporative (blue 
color) and midvelocity (red color) multiplicities of p, a 
and IMFs in the reaction 93 Nb + 93 Nb at 38 AMcV (the 
three cases correspond to dominant evaporation, com- 
parable components and dominant midvclocity emission, 
respectively); the black curves show the corresponding 
total multiplicities. The reference (shown by full dots) 
is always the analysis with the choices adopted in this 
paper, namely evaporation in the angular range 0°-45°, 
extrapolation - based on Monte Carlo simulations - for 
non-zero spin source and applied corrections for recoil ef- 
fects. Variations are then made one at a time: either a 
zero-spin source is assumed (triangles), or the recoil cor- 
rections are switched off (squares), or the angular range 
is reduced to 0°~30° (stars) or even enlarged to 0°-90° 
(open dots). 

Being the decomposition between M eV ap and M m idv 
obtained with a subtraction procedure, these different 
choices induce larger relative variations in the weaker 
component: thus, in the evaporative component they are 
largest for the IMFs (which arc predominantly emitted 
at midvelocity) and, in the midvelocity component, they 
are rather large for protons at 38 A MeV and even larger 
- this time for all particles - at the lowest bombarding 
energies (where evaporation dominates). 

In general one can observe that: 

• The more conservative angular range 0°-30° gives 
values of A4 C vap (-Mmidv) which are just slightly 
smaller (larger) than those obtained with 0°-45°. 
Only for IMFs there is a sizable reduction (about 
-30%) of the already quite small A^evap component 
(while the corresponding increase of M m idv is neg- 
ligible). 
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FIG. 17: (color) Total (black), evaporative (blue) and mid- 
velocity (red) multiplicities for p, a and IMFs (upper, middle 
and lower panels, respectively) in 93 Nb+ 93 Nb at 38 A MeV. 
The results of the analysis adopted in the paper (full dots) are 
compared with those obtained from slightly different analysis 
procedures (see text). 



• With the commonly used procedure of estimat- 
ing the evaporation by taking twice the particles 
emitted in 0°-90°, one systematically overestimates 
■Mevap (and hence underestimates M m idv), for all 
particles and at all bombarding energies. For pro- 
tons, the variation is negligible on M cva , p and small 
but sizable on A4 m idv (especially at high TKEL); 
for a the effect is of the order of 30-40% and com- 
parable on both components; for IMFs the effect 
is huge (about a factor of 2-3) on the weak A'fevap 
component and much smaller - but still sizable, 
about 20-30% - on the dominating A^midv In- 
deed, for IMFs this procedure gives the largest of 
all considered variations. 

• Extrapolating the angular distribution of evapo- 
rated particles with a cos 9 plf shape generally 
overestimates A'fevap (and underestimates A^midv)- 
However, at small TKEL (where the spin is likely 
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to be small), it produces really negligible variations 
in all cases. At large TKEL, there is a moderate 
effect (of the order of 10-20%) only on the weaker 
component. 

• As explained in Sec. IIII El applying no recoil cor- 
rection causes an overestimation of TKEL, namely 
a given bin erroneously includes some events (which 
should actually be classified in a bin at lower 
TKEL) and misses others (which get classified in 
another bin at larger TKEL). The net effect of 
this gain and loss is determined by the TKEL de- 
pendence of the experimental particle yield, i.e., 
yield of measured events times event multiplicity 
for that particle). The rise of all multiplicities is 
steepest at small TKEL and tends to flatten at 
large TKEL, while the measured event yield doesn't 
change more than a factor of two over the whole 



TKEL range. Therefore the effect of the recoil cor- 
rection is largest at small TKEL, where it deter- 
mines an increase of A^evap - of course larger for 
heavier particles - from 10-20% for protons to a 
factor of about 2 and 3 for a particles and IMFs, 
respectively (the corresponding decrease of M m idv 
is about 70% for protons, a factor of 2 for a and 
negligible for IMFs). 

In the Nb+Nb system at 38 A MeV, the effects of the 
considered analysis variations are most visible in the 
PLF-multiplicities of IMFs at large beam energies, and in 
the midvelocity multiplicities of LCPs at low beam ener- 
gies and large TKEL. For all investigated reactions, these 
effects are schematically represented by the error bars in 
Figs. and and they should be kept in mind when 
comparing the results of different experiments obtained 
with different analysis procedures. 
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